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We analyze the meson correlator in the spatial direction at finite temperature. To achieve fine resolution in 
the spatial direction, we use an anisotropic lattice with the standard Wilson plaquette gauge action and the 0(a) 
improved Wilson quark action. Below and above T c , properties of correlators are investigated by two methods: 
fits with ansatz for the spectral function, and direct reconstruction of the spectral function using the maximum 
entropy method. 



1. Introduction 

The extraction of the spectral functions from 
the correlators measured in lattice QCD simu- 
lation is of great importance to understand the 
QCD phase transition at finite temperature, both 
in theoretical and phenomenological points of 
view. Our goal is to investigate the both of spa- 
tial and temporal structure of hadron correlators 
near the deconfining transition Jl[. Particularly 
for the latter, one needs to develop techniques to 
extract detailed and reliable information from the 
correlators. On the other hand, the correlators in 
the spatial direction are rather easy to measure 
precisely, and therefore are expected to serve as a 
good test ground for the method. In this report, 
we focus on the spatial meson correlators at finite 
temperature, and develop procedures to analyze 
them. These analyses will give us insights also for 
the static screening structure of QCD at T > 0, 
which still remains not well understood ||. 

We introduce a spectral function for the corre- 
lator propagating in the spatial direction in anal- 
ogous to the ordinary spectral function obtained 
from the correlator propagating in the temporal 
direction. In order to reconstruct this spectral 
function from the lattice data, we apply two kinds 
of analysis procedures: (a). Fit with ansatz on 
the shape of spectral function. At present stage, 
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we assume that the spectral function is repre- 
sented by a sum of several strong peaks, i. e. delta 
functions, and fit the data to this form to obtain 
the screening masses corresponding to the peaks. 
We adopt single, two, and three pole(s) fittings, 
(b). Reconstruction of spectral function from 
the correlator by the maximum entropy method 
(MEM) H. Comparing the results of these two 
methods, we discuss their efficiency and reliabil- 
ity. 

For these analysis, it is better to have as many 
points of correlators as possible. We use the 
anisotropic lattice which has small lattice spacing 
in z-direction to measure the spatial hadron cor- 
relators in high resolution. Although one needs 
additional effort for the calibration of anisotropy 
parameters, once it has been done, the system- 
atic uncertainties due to the anisotropy can be 
controlled within reasonable accuracy B. In the 
following, we show only the results for the pseu- 
doscalar meson correlators, while the procedures 
are applicable for other channel in the same man- 
ner. 

2. Numerical Simulations and Results 

Our numerical simulations are carried out on 
quenched lattices of sizes 12 2 x 96 x N t where 
N t =12, 5, 3 and 2 which roughly correspond 
to the temperatures ~ 0, 0.8T C , 1.3T C and 2T C , 
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respectively. The gauge configurations are gener- 
ated with the Wilson plaquette action at j3 = 5.75 
with the anisotropy £ = a/a z — 4, where a z 
and a are the lattice spacings in z-direction and 
other three directions Q|. With these parame- 
ters, the lattice cutoff set by the string tension is 
a^ 1 ~ 4 GeV. Configuration numbers are 400 at 
T ~ and 510 at other T's. The quark action 
is the 0(a) improved Wilson action Q, at three 
values of hopping parameter, k =0.124, 0.122 
and 0.120, which cover the range of quark mass 
m q — (0.5 ~ 1.5)m s . The bare quark anisotropy 
is set as 7f =, which is determined using the 
meson dispersion relation from the lattice Klein- 
Gordon type action. 

In the following, we show the results for the 
pseudoscalar meson channel at k = 0.120 (m q ~ 
1.5m s ). The fit is performed with the fitting re- 
gion [L min , L max ], where L max is fixed to 48 and 
L m in is varied. The screening masses as the re- 
sults of single, double and triple pole fits are dis- 
played in the figures for each value of L min . Ex- 
istence of a region where the masses are stable 
indicates that the assumed form of fit is a good 
representation of the spectral function. On the 
other hand, absence of such a region signals that 
the spectral function is not approximated by such 
a form. In both cases, we can extract the spectral 
function reconstructed by MEM as an alternative 
procedure. In MEM analysis we use the default 
model function mo(u) = 16w 2 , frequency region 

< uj < 2.0 and the fitting region of correlators 

1 < t < 40. 

Let us start with the zero temperature. The 
upper panel of Figure [l] shows the results fitting 
analysis. For each fit (i.e. number of poles) , there 
is a region where the masses do not change within 
errors. For the ground and first excited states, 
three fits give consistent masses. This shows that 
the assumed form of spectral function (the sum 
of poles) well represent the correlator, in accord 
with the physical expectation. The spectral func- 
tion reconstructed by MEM also exhibits consis- 
tent position of peaks, as shown in the lower panel 
of Fig. 0. Since the masses from the fits shows 
clear plateaus, and observing the errors associ- 
ated to the regions of spectral function, the ob- 
served widths of peaks may not be practical ones 
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Figure 1. Results at T = 0. The upper panel 
is the result of fits. The lower panel shows the 
spectral function obtained using MEM. 



but artifacts through the statistical and system- 
atic errors MEM. 

Now we apply the same procedures to the cor- 
relators at finite temperature. The result of fit 
method at 0.8T C is almost the same as at T = 0. 
Therefore there seems no significant change in 
the spectral functions at T — 0.8T C from that at 
T = 0. On the other hand, as shown in Figure 
the result of MEM shows a deference: there is two 
peaks at the position of the first excited state ob- 
tained by the fits. This is unphysical structure 
can be understood by the failure of MEM in re- 
constructing the spectral function, probably due 
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Figure 2. Result for the spectral function at T = 
0.8T C . 



to lack of statistics. In order to detect the struc- 
ture of spectral function in high frequency re- 
gion from MEM, one would probably need larger 
statistics that the present one which of O(500). 

Figure || shows the results at 1.3T C . We again 
find the regions of constant masses for each fit and 
a consistency in the masses of the ground and first 
excited states. Result of spectral function also 
have peaks at consistent positions. Thus we find 
that the sum of poles are a good approximation 
of the spectral function. The ground state mass 
is slightly less than the twice lowest Matsubara 
frequency. For the quantitative analysis, however, 
present lattice is so coarse that the number of 
temporal degrees of freedom is just three. We 
also observe similar result at T = 2T C . 

We conclude that these two procedures are 
promising as useful ways to extract the structure 
of the spectral function, although one need care- 
ful treatment particularly for MEM. Fits assum- 
ing various forms for the spectral function serve 
complementary way of analysis, and supply the 
reliability for the result of MEM. Such a study 
is particularly important for the analysis of tem- 
poral correlators, since the degrees of freedom is 
restricted even with anisotropic lattices and the 
data are inevitably in the high frequency region. 
In this case, since the spectral function have phys- 
ical width, one also need to try other fitting forms, 
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Figure 3. Results at T = 1.3T C . 

such as the peaks with widths. To obtain more 
quantitative result on the screening mass and the 
spectral function for the spatial correlators, we 
should use smaller lattice spacing as well as higher 
statistics. 
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